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Abstract. NLSP - LSP decays could have dramatic influence on SUSY phenomenology at LHC. 
NLSP could have significant lifetime and could be charged. In at least two scenarios detectors must 
be used in a special way. They were not optimized for detection of heavy (semi)stable charged 
particles and decaying in flight (neutral or charged) NLSPs. During the last decade both ATLAS 
and CMS collaboration have developed strategies which allow for effective search within such 
scenarios. 

PACS. 14.80.Ly Supersymmetric partners of known particles - 29.00.00 Experimental methods 
and instrumentation for elementary-particle and nuclear physics 



1 Introduction 

Within Gauge-Mediated Supersymmetry Breaking 
models gravitino is the Lightest Supersymmetry Part- 
ner (LSP) whereas neutralino or stau plays the role of 
the Next to Lightest Supersymmetry Partner (NLSP). 
NLSP decays to its Standard Model partner and grav- 
itino with a lifetime depending on the scale of the 
SUSY breaking. Detection of the NLSP and determi- 
nation of its properties, in particular lifetime, can be 
of crucial importance for the physics program at LHC. 

Minimal GMSB model is defined by six parameters: 
A - effective SUSY mass scale, N - number of mes- 
senger generations, M - messenger mass scale, tan/3, 
sgn/Lt and Cg - ratio of the intrinsic SUSY breaking 
scale to messenger SB scale, governing goldstino cou- 
pling and hence - NLSP lifetime. 



2 Long-lived charged NLSP at CMS 

For large N, typically right sleptons are the lightest su- 
persymmetric partners of SM particles. If, in addition, 
tan f3 has not too low value, fi is the NLSP which could 
have arbitrary long lifetime. In this way one obtains 
benchmark scenario in which in every supersymmetric 
event cascade decays ends on pair of stable massive 
charged and not strongly interacting particles. These 
particles propagate through the detector like muons 
but with velocity /3 smaller than 1. Because of that 
their specific ionization is greater than for MIP and 
they arrive at given detector layer with a time de- 
lay. Both facts could be used to distinguish them from 
muons. 

There are other theoretical scenarios with heavy 
stable charged particles. ATLAS and CMS searches for 



them were reported by Shikma Bressler in the same 
session [1]. I will cover only some details concerning 
search for staus using TOF method in the CMS de- 
tector. The first full detector simulation (based on 
GEANT3) analysis were performed already a decade 
ago [2l[3]. However specific ionization were not cor- 
rectly simulated then. Because of that the analysis was 
based on originally developed method of TOF mea- 
surement by drift tubes of the barrel muon system of 
the CMS. 

Here we present an update [4] of that analysis per- 
formed using OSCAR-ORCA version of CMS detector 
simulation (based on GEANT4), the same as used for 
CMS Physics TDR 0. 



2.1 The TOF method 

The barrel muon system of the CMS detector consist 
of four concentric muon stations inserted in the return 
yoke of the CMS solenoid. In each station there are 
three super-layers (SL) of four layers of drift tubes 
(DT) each. Two of SL measure R(j) coordinate and 
one z coordinate (there is no z SL in the outermost 
station). The tubes in each SL are staggered by half 
a tube. An average track pass alternatively on left and 
right side of the sensitive wire. 

In a given SL hits due to muon should align if tim- 
ing is correct whereas hits due to delayed particle do 
not align, they are shifted backward from the wire by 
dx and form a zig-zag pattern. For each hit 
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Table 1. Number of events at different stages of the se- 
lection for 1/fb. Columns correspond to: presel. - trigger 
and > 80 GeV; quality - additional requirements at the 
pattern recognition stage that significantly reduce tails of 
f3~^ estimate (see ref. [4] for details); select. - requirements 
described in the text; - highlighted area in the Fig- 
ure [1] 
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Fig. 1. Scatter plot /3 ^ versus momentum 



The exact formula used in the analysis to obtain 
estimate from the zig-zag pattern is described else- 
where [4] but it could be simplified as follows: 



P I'drift N ^ Li 

where Li is the flight distance, Wdrift is the drift veloc- 
ity, SI = \xf^ - xY"°\ - 1x1°""° - and is the 
local track element position at layer i. 



2.2 Results 

Two points from the SPS7 line 6 : yl = 50 TeV and 
yl = 80 TeV were chosen for the full detector simula- 
tion. The mass of fi particle is equal to 152.31 GeV 
and 242.93 GeV respectively and the cross section is 
Ipb and 0.1 pb. In the analysis a high mass (Af^^ > 
llOGey) pair of energetic muons {pt^ > 60 GeV) were 
searched for. Since most of the signal events contains 
products of a cascade decays of squarks or gluinos a 
high effective mass of the event were required {Mos > 
360 GeV) as well. Taking this into account the fol- 
lowing background sources of muon pairs were con- 
sidered: Drell-Yan above resonance, tt and double 
vector boson production. Events were triggered on sin- 
gle muons with a default threshold of {pt^ > 80 GeV). 

Breakdown of number of events at different stages 
of the selection is given in the Table [T] A scatter plot 
versus momentum for heavier mass point and C = 
4/fb is shown in the Figure [1] where clear separation 
of signal - staus and background - muons could be 
seen. The last column in the Table [T] correspond to 
highlighted area in that Figure. 

The upper limit for number of expected background 
events in the region above the highest horizontal line 
in the figure (/3~^ > 1 + 3ap~i = 1.132) was evaluated 
to be 0.05 events for 1/fb which means that the 5a 
discovery could be claimed with 8 signal events. This 
correspond to integrated luminosity needed for a dis- 
covery £ = 52/pb for yl = 50 TeV and £ = 667/pb for 
A = 80 TeV. 



It should be stressed, however, that real perfor- 
mance of the analysis could not be evaluated without 
checking the method against muons from real data. 

Another question is how well the fi mass could be 
determined. To test this 1000 pseudo-experiments cor- 
responding to integrated luminosity £ = 0.5/fb for the 
lower mass point {£ = 4/fb for the higher mass point) 
were performed. The result are the following 
M|^*- = {153.2 ± 1.6(stat.) ± 0.9(syst.)}GeV 

for generated mass M5°"' — 152.31 GeV and 

= {243.2 ± 3.2(stat.) ± 1.4(syst.)}GeV 

for generated mass AfS°"' = 242.93 GeV. The test con- 

firms that the method is well suited for the search for 
stable staus. 

3 Non-pointing photons in ATLAS 

Within GMSB if N is small typically x? is the NLSP. 
Its lifetime is free parameter. By measuring lifetime 
it is possible to estimate the scale of supersymmetry 
breaking which is otherwise not accessible experimen- 
tally. 

In the ATLAS analysis [7] very interesting tech- 
nique was developed which allows to determine the 
masses of the slepton and neutralino from events with 
a lepton and converted photon arising from the cas- 
cade decay £ £xi — > £jG. 

The topology of Xi decay inside ATLAS detector if 
the decay length is of the order of 1 m is shown in the 
Figure [H In the case of photon conversion the angle a 
could be precisely measured. ATLAS EMCAL allows 
also for precise determination of photon arrival time 
t.y. Since the distance L = \OA\ is also known the 
angle ip could be determined [7] 

1 — , ^ ct^ — L cos a 

cos tp = — where t, = — . 

1 -t- 4"^ _Lsma 

If both angles a and -0 are know for at least two 
reconstructed £ Ixi ~^ cascade decays it is 
possible to determine both £ and Xi masses. 
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Fig. 4. Average {to /"fx) {^) number of lepton non- 
Fig. 2. The topology of x? decay inside ATLAS detector. pointing photon pairs (b). 
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Fig. 3. Reliability of determination (a) and to /"fx dis- 
tributions (b). 



Knowing these masses it is possible to analytically 
calculate decay time and momentum of neutralino us- 
ing the ECAL data and lepton momentum only, for 
each lepton non-pointing (converted or not) photon 
pair from slepton decay. 

This idea was tested with fast simulation of ATLAS 
detector for GMSB point: A = 90TeV, M = 500 TeV, 

= 1, tan/3 5, /i > 0. Neutralino and slepton 
masses for this point are Af(x°) = 117 GeV, M{£ji) = 
162 GeV 

Longitudinal segmentation of the EMCAL of AT- 
LAS allows for precise determination of the polar angle 
of a non-pointing photon with very good resolution of 
O.OQ/ y/E^ /GcV. More over, the arrival time could be 
measured with 100 ps resolution. This in turn allows 
for determination of Xi decay time t^i. In the Figure [3] 
a correlation between reconstructed and generated to 
is shown as well as distributions of i^) /j^ for three val- 
ues of CT. In the Figure[4]the average {to Mx) ^'^'^ num- 
ber of lepton non-pointing photon pairs in function of 
generated cr is shown after the following selection: 

> 30 GeV, A0 > 0.2, At^ > Ins M^s > 400 GeV, 
£'^""" > O.lMcff for integrated luminosity C = 13.9/fb. 

Both values shown in the Figure|3]could be parame- 
terized and use for Xi ct determination. The reliability 
of both methods and absolute resolution (error bars) 
in function of generated cr is shown in the Figure [5] 



The relative resolution of cr ranges from 3% to 17% 
for average proper time method {tn/lx) ^'^"^ form 3% 
to 6% for lepton non-pointing 7 pair counting method. 
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Fig. 5. Reliability of cr determination methods. 

4 Non-pointing photons in CMS 

A feasibility study of x? hfetime determination using 
CMS detector was done long ago [3]. However, only 
recently and only part of the original proposal was 
transformed into full detector simulation analysis [4]. 
The method is based on a difference between shapes 
of energy distributions among CMS ECAL crystals for 
pointing and non-pointing photons (see Fig. [6]) . 

These energy distributions could be characterized 
by the covariance matrix. Variances along major and 
minor axes are given by 



»S' major 



J(500-5'^„)2+4^; 
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A measure of the elongation of the deposit is an asym- 
metry 
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Another useful variable is the angle a (see Fig. [6]) be- 
tween the major axis and the cj) axis which could be 
used to eliminate background due to converted pho- 
tons for which a « and geometrical bias which gives 
an excess of deposits with large A for pointing photons 
for a « 0,±f . 

For the full simulation a point from SPSS line [B] 
with A — 140 TeV was chosen and 6 different x^ values 
of cr were generated. The x^ mass is 192.69 GeV and 
the cross section is 0.5 pb for this point. 

The analysis is based on a search for energetic pho- 
tons which are accompanied by at least four hard jets 
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Fig. 6. Energy distribution in tlie ECAL crystals for point- 
ing (left) and non-pointing (riglit) piioton. 

Table 2. Number of events at different stages of the se- 
lection for signal and background for jC — 10/ fb. Prese- 
lection: single 7 trigger, > 80 GeV; event selection: 
< 1.7, Zi0(jj) > 20°, pj? > 50 GeV, p^''' > 160 GeV 
non-pointing 7 selection: A > ^, a ^ 0, ^ pointing 



7 selection: NOT non-pointing, 
220 GeV. 



> 



60 GeV, p^""" 



> 



dataset 


presel. 


select. 


non-p. 


point. 


CT = 


2104.28 


402.98 


2.96 


289.39 


CT — 25 cm 


2061.69 


379.71 


29.10 


243.76 


cr — 50 cm 


1948.33 


361.60 


45.87 


215.88 


CT — 100 cm 


1564.12 


298.12 


56.03 


182.99 


CT = 200 cm 


1037.60 


166.66 


34.64 


99.39 


cr — 400 cm 


645.77 


114.89 


26.16 


67.19 


rZjets 


6137.93 


0.65 


0.00 


0.37 


rWjets 


8301.93 


2.76 


0.00 


1.46 


EQCD 


787797.19 


54.90 


0.27 


2.32 


VVjetsJncl. 


836.35 


O.QQ 


0.00 


0.00 


ttbarancl. 


4662.89 


16.34 


0.00 


6.13 


rtotal 


807736.31 


74.65 


0.27 


10.27 
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Fig. 7. £ for 5(t discov- 
ery versus cr. 



Fig. 8. Estimated cr 
versus generated cr. 



and missing energy, a consequence of cascade decays 
of squarks and gluinos. A breakdown of the number 
of events at different stages of the selection for signal 
and background samples is given in the Table [2l 

The fraction of selected photons classified as non- 
pointing improves with growing cr but at the same 
time more and more neutralinos decays outside ECAL. 
Non-pointing selection is almost background free. 

An estimate of the minimal amount of integrated 
luminosity needed for Scr discovery versus simulated 
CT is plotted in the Figure [T] 

An integrated luminosity C — 2/fb is sufficient to 
claim discovery for cr < 100 cm and C — 3.5/fb for 
cr = 400 cm. 



To test CMS ECAL capability to estimate Xi life- 
time 100 likelihood fits of neighborhood cr and back- 
ground to the given cr + background was done. The 
results are shown in the Figure[8l The precision of such 
procedure range from 15% to 40%. It will improve if 
more cr points will be used, but one should remem- 
ber that it depends on the knowledge of the shapes of 
distributions of variables used in the fit. These shapes 
depend (weakly) on other than cr parameters of the 
model. 



5 Conclusions 

Gauge-Mediated Supersymmetry Breaking model is a 
generic framework for long-lived, heavy, charged, not 
strongly interacting particles and decaying in flight 
heavy (charged or neutral) particles. 

During last decade both Atlas and CMS have in- 
vented methods to tackle with such signatures, despite 
the fact that detectors had not been designed for them. 

For long lived charged particles both detectors use 
synergy of specific ionization measurement and TOF 
method (for the full picture see also Ref. TJ). The ulti- 
mate goal before LHC startup should be development 
of model independent methods which could be tested 
on cosmics, Z f-^ f-^ candle and other sources of 
energetic muons. 

For non-pointing photons signature ATLAS and 
CMS successfully use their electromagnetic calorime- 
ters. However, the ultimate resolution for lifetime de- 
termination of decaying in fiight particles could be 
obtained only with tracking systems. Performance of 
such methods require full detector simulation with full 
background turned on. 

There is just enough time to do that before LHC 
startup. 
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